APPENDIX 1 



%%% 

function B = Initcluster(X,m) 

% Get an initial base codebook B at random from input data 
% INPUTS 

% X = input data: each column is a RGB 'vector' 
% m = number of base codevectors 

% OUTPUTS 

% B = base codevector matix 

%need to duplicate some columns of B if X is small 

[n,N] = size(X); 

if(N>m) 

replace = 0; 
else 

replace = 1; 
end 

%track what inputs put in B, so no duplication for X large 
chosen = zeros(l ,N); 
B = zeros(n,m); 
for i=l :m 
draw = floor(N*rand +1); 
if(~replace) 
while(chosen(draw)) 

draw = floor(N*rand +1); 
end 
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end 

B(:,i) = X(:,draw); 
chosen(draw) = 1 ; 
end 

%%% 

function [B,Bq,d,iters,rd,nX] = lgbo(X,m,base) 
% INPUTS 

% X = input data: each column is a RGB 'vector' 

% m = number of codevectors (columns in C) 

% base = number of base codevectors (columns of B) 

% OUTPUTS 

% B = base codevector matrix 

% Bq = quantized base codevector matrix 

% d = distortion of X when replaced with chosen codevectors 

% iters = # of iterations to reach convergence 

[d,N] = size(X); % d = dimension, N = blocksize 

%random initialization 

Init = initcluster(X,base); % choose random initial set 

if base==2 %1D (linear interpolation) 

W = [l:(-l/(m-l)):0; 0:(l/(m-l)):l]; %weight matrix 
elseif base=3 %2D 

W=[l,0,0;0,l,0;.25,.25,.5;-.25,-.25,1.5; ... 
1,-1,1;-1,1,1;0,-1,2;-1,0 5 2;]'; 
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elseif base=4 %3D 

W=[eye(4),. 25 * [2,2, 1,-1;2,- 1,2,1; 1,2,- 1,2;- 1,1, 2,2]]; 
%W=[eye(4),.125*[3,2,2,l;2,3,2,l;3,l,2,2;l,2,3,2]]; 
end 

%variable initialization 
stoppingeps = l.e-5; 
vi = ones(l,m); 

index=zeros(l ,N);mind=index; 
cumdist = Inf; lastdist = 0; 

C=Init*W; %interpolate codevectors 

# iterate until convergence 
while(abs(cumdist - lastdist) > stoppingeps) 

lastdist = cumdist; 

cumdist = 0; 

iters=iters+l; 

while(l) %1 iteration that repeats if B goes singular 
% step (A) 

% form Voronoi regions: for each input, 

% determine which centroid it is closest to 
fori=l:N 
V=X(:,i(vi))-C; 

nm=sum(V.*V); %Euclidean distance squared (MSE) 
%input i's closest codevector 
[mind(i),index(i)]=min(nm) ; 
end 

cumdist = sum(mind); 



28 



%find diagonal matrix N 
for j=l:m 

n(j) = sum(index=j); 

5 end 

% check if B is singular 
% and force it to be non-singular 
num==sum(n~=0) ; 
10 ifnum<base %fewer than base non-zero! 
nz = fmd(n); 

if (base-num)==l %base==2 
[jy,ji]=max(niind); 
Init(:,l :(base-l))=B(:,nz); 
15 Init(:,base)=X(:,ji); 
else 

D y J i]=sort(mind) ; 
Init(:,l :num)=B(:,nz); 
jl=N; 

20 forjk=(num+l):base 
mit(:,jk)=X(:,ji(jl)); 

while(sum(abs(X(:,jiGl))-X(:jiai-l))))=0) 

jHl-1; 
end 

25 jl=jl-l; 

end% forjk 
end % else 

C=hiit*W; %start all over if hit a singular matrix 

else 

30 break; 
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end % if not singular then end while loop 
end%while(l) 

s = zeros(d,m); %get sum matrix S 
for j=l:m 

ifnG)=l 

s(: j) = X(:,index=j); 
elseif n(j) 

s( : j )=sum(X(: ,index==j )')' ; 
end 

end %for j=l:m 

Init = s*W'*inv(W*diag(n)*W'); %new base codevector matrix B 
C=Init*W; %new codevector matrix C 

end %% while not converged 

nX = B(:,index); % save re-constructed block 
d = cumdist; % total block distortion 

end % function 
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